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Abstract 


Pileus clouds form where humid, stably stratified air is mechanically displaced 
vertically ahead of rising convection. This paper describes convective forma- 
tion of pileus cloud in the tropopause transition layer (TTL), and explores a 
possible link to the formation of long-lasting cirrus at cold temperatures. In-situ 
measurements from off the coast of Honduras during the July 2002 CRYSTAL- 
FACE experiment show an example of TTL cirrus associated with, and pene- 
trated by, deep convection. The cirrus was enriched with total water compared 
to its surroundings, but composed of extremely small ice crystals with effective 
radii between 2 and 4 pm. Through gravity wave analysis, and intercomparison 
of measured and simulated cloud microphysics, it is argued that the TTL cirrus 
in this case originated neither from convectively-forced gravity wave motions 
nor environmental mixing alone. Rather, it is hypothesized that some combina- 
tion was involved in which, first, convection forced pileus cloud to form from 
TTL air; second, it punctured the pileus layer, contributing larger ice crystals 
through interfacial mixing; third, the addition of condensate inhibited evapora- 
tion of the original pileus ice crystals in the warm phase of the ensuing gravity 
wave; fourth, through successive pulses, deep convection formed the observed 
layer of TTL cirrus. While the general incidence and longevity of pileus cloud 
remains unknown, in-situ measurements, and satellite-based Microwave Limb 
Sounder retrievals, suggest that much of the tropical TTL is sufficiently hu- 
mid to be susceptible to its formation. Where these clouds form and persist, 
there is potential for an irreversible repartition from water vapor to ice at cold 
temperatures. 
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1 Introduction 


At low latitudes, the tropopause transition layer (TTL) represents the layer between ap- 
proximately 14 to 18 km altitude through which tropospheric air enters the stratosphere 
(Rosenlof, 2003). It is characterized by high vertical gradients in water concentration and 
a local minimum in temperature. The mechanisms governing distributions of water vapor 
30 and ice in the TTL have attracted considerable interest. This is primarily because these 
species contribute to the planetary greenhouse mechanism by strongly absorbing terrestrial 
radiation at temperatures typically 100 K colder than the surface. 

This paper attempts to describe a process in the TTL that has been largely overlooked, 
but is visibly associated with a repartition of water from vapor to ice. When clear, moist, 
35 stratified air is pushed upward above rising deep convection, it cools adiabatically. The 
activity of solution aerosols in the air increases, sometimes to the point that they freeze 
homogeneously. The newly formed ice crystals subsequently grow by the diffusion of water 
vapor to their surface. The process is visible as a tenuous veil that forms over the tops of 
convective cloud turrets. These clouds have been called pileus, from the Latin word for the 
40 felt caps that were worn by freed slaves. 

The motivation for studying cloud formation of any type in the TTL is that the effects of 
water molecules on atmospheric processes is different in the condensed phase than it is in 
the gas phase. For one, condensate settles, and can dessicate air. Second, when molecular 
dipoles, such as H 2 O, arc placed in close proximity, as they are in condensate, they respond 
45 to electromagnetic radiation by vibrating in phase, greatly amplifying the interaction with 
electromagnetic radiation. 

In the case of pileus clouds, however, normally it would be assumed that their forma- 
tion is associated only with a transient influence to water vapor concentrations and radiative 
fluxes - once the pressure forcing is gone, air returns to its previously clear initial atmo- 
50 spheric state. Flere, we examine this assumption more closely using photography, in situ 
measurements, and numerical simulations. The results suggest that, where the TTL is ini- 
tially supersaturated, or pileus is punctured and mixes with deep convection, pileus cloud 
may linger. If so, an irreversible repart ion from vapor to condensate occurs within TTL air. 
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2 Photography 


55 For illustration, we show two examples of pileus clouds forming near the tropopause. The 
first shows pileus forming in ’’stacked” layers above a deep convective storm that developed 
over land in the Tropical Western Pacific region near Darwin, Australia (12° S, 130° E) 
(Fig. Al). At the Darwin airport, fifteen minutes after this photograph was taken, there was 
extensive lightning, an unusually intense surface gust front of 83 km hr 1 , a temperature 
60 drop of 10°C, and the precipitation intensity reached 80 mm hr -1 . It is impossible to say 
with certainty, but the vigor of the storm is suggestive that, near the time the photograph 
was taken, the convective tops were near the tropopause, which local soundings indicate 
was located at 17.8 km and 186 K. According to the apparent height of the convective cloud 
in the photograph, the convection had pushed the moist layer that formed the pileus cloud 
65 about 10% (2 km) higher. 

A similar example of pileus formation above deep convection has been described by 
Scorer (1972), based on observations in Kenya. A series of four photographs showed deep 
convection forming then puncturing pileus cloud. Once the convection subsided, the pileus 
lingered and spread. At the same time, anvil outflow from the deep convection spread as 
70 a separate cloud layer beneath the newly formed cirrus. Photographs in Fig. A2, which 
were taken over Fouisiana at an altitude of 12 km, appeal - to show a similar process. Tur- 
bulent cirrus anvil outflow spreading from deep convection was apparently accompanied by 
two thinner, laminar cloud layers at higher altitudes. It is possible the highest layer was 
pre-existing. However, its horizontal extent, and the presence of wave-like features with 
75 dimensions similar to the width of the deep convective turret, suggests the layer originated 
as pileus cloud forced by convective uplift. The lower layer, which was more clearly pileus, 
enveloped the convective turret in a characteristic skullcap veil. Six minutes later, this pileus 
layer had been punctured by the turret, and spread laterally above the anvil. Unfortunately, 
no photographs were taken at later stages of development when the convection had sub- 
80 sided. However, at least at its edges, this lower layer appeal's to have survived to the point 
where it was evolving as TTF cirrus independent from the convection that forced it. 

3 Formation of TTL cirrus near deep convection 

The photographs show that convection can induce cirrus formation at cold temperatures and 
low latitudes; observations appeal" to indicate that sometimes overshooting deep convection 
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85 can puncture a pileus cloud it forms. This is notable, because interfacial turbulent mixing 
is generated by shear between convective cloud and its environment (Grabowski and Clark, 
1991, 1993), and humidity becomes enhanced in clear air up to several cloud radii distant 
(Perry and Hobbs, 1996; Lu et ab, 2002). We now discuss whether a possible exchange 
of air between convection and its surroundings may play a role in the evolution of a pileus 
90 cloud. 

While we know of no in situ observations in the TTL that were conclusively made within 
pileus cloud, we describe here measurements of TTL cirrus obtained from the NASA WB- 
57F aircraft during the July 2002 CRYSTAL-FACE campaign over Southern Florida and 
near Honduras. During CRYSTAL-FACE, measurements of ice concentrations were ob- 
95 tained with a Cloud and Aerosol Particle Spectrometer (CAPS) (Baumgardner et al., 2002), 
water vapor w and total water concentrations Q with the Harvard Water Probe (Weinstock 
et al., 1994), the optical extinction coefficient f3 from a Cloud Integrating Nephelometer 
(CIN) (Gerber et al., 2000), total water isotopic concentrations HDO and H 2 O with the 
Aircraft Laser Infrared Absorption Spectrometer (ALIAS) (Webster et al., 1994), and tem- 
1 00 perature, pressure, and wind speed from the Meteorological Measurement System (MMS) 
(Scott et al., 1990). 

From observations on six flights, Garrett et al. (2004) noted that cirrus had formed in the 
TTL directly above, but separate from anvil outflow from deep convection. In these cases, 
the anvil and TTL cirrus did not appear to have formed independently. Cloud boundaries 
1 05 established from aircraft and satellite indicated that the TTL cirrus had nearly the same hori- 
zontal extent and location as the anvil. Also, mixing ratios of nitrogen oxide (NO), potential 
temperature ( 0 ), and total water in the TTL cirrus were intermediate to those observed in 
clear TTL air and those in anvil cirrus. Thus, it was interpreted that the formation of TTL 
cirrus involved a fractional contribution of air from deep convection. This contribution 
1 1 0 ranged up to 0.5. In one TTL cirrus layer suitable for such analysis, measured on 21 July, 
spectral analysis of the temperature field showed that it contained a distinct monochromatic 
gravity wave signature not seen in surrounding clear air. The wavelength of this feature was 
~2 km - similar to the dimension of deep convective turrets - and its ampltitude was several 
hundred meters. The apparent combination of mixing and monochromatic wave motion led 
115 to the speculation that the TTL cirrus was forced by deep convection, but began as pileus 
cloud that had mixed with the convective turret. 

The extent of mixing between different airmasses can be inferred from measurements of 


4 



such conserved variables as enthalpy and total water. In addition, the airmass origins can 
be inferred from relative concentrations of isotopes of total water, since these fractionate 
1 20 according to their condensation temperature - a proxy for their altitude (Kuang et ah, 2003; 
Webster and Heymsfield, 2003), and in total water, they mix linearly. According to Rayleigh 
distillation, atmospheric water becomes progressively lighter as an air parcel is cooled; 
heavier isotopes preferentially condense and precipitate. In the absence of any mixing, the 
depletion of HDO relative to H 2 O (S— HDO) would range from -8.6% above the ocean, to 
1 25 ~-95 % at the coldest tropical tropopause. Recent measurements have shown average values 
of 5 — HDO in the TTL arc in fact near -64%, which argues that the origins of TTL moisture 
are determined by condensation processes well below the tropopause (Kuang et ah, 2003). 
However, Webster and Heymsfield (2003) showed that this value is only an average, and 
TTL values of 6 — HDO arc actually highly variable, ranging from -90 to ~0% at horizontal 
1 30 scales down to ~3 km. 

An illustration of a role for mixing in TTL cirrus formation is shown in Figure A3, from 
a vertical descent during CRYSTAL-FACE by the WB-57F, first through TTL cirrus, and 
then followed by anvil cirrus. The aircraft flew through what, from satellite imagery, looked 
like a continuous anvil cloud shield, entering the TTL cirrus about 40 km downwind from 
135 where the anvil leading edge was located, and entering the anvil layer a further 20 km 
downwind. Values of 4- HDO in the TTL cirrus were lower than those observed in anvil 
air beneath it, but higher than surrounding TTL air, which suggests the TTL cirrus was 
formed by mixing convective air into the TTL. It might be argued that the anvil cirrus and 
TTL cirrus were equally derived from deep convection: i.e. that the TTL cirrus formed 
1 40 simply from detrainment of the same convective air mass that formed the anvil; HDO was 
lower in the TTL cirrus simply because it had been depleted by precipitation. However, the 
potential temperature in the TTL cirrus was about 3 K colder than air immediately above 
and below, and about 10 K warmer than the anvil air. Also, total water concentrations in the 
TTL cirrus were intermediate to those seen in surrounding TTL air and anvil cirrus below. 
1 45 Because enthalpy, 6— HDO, and total water mix linearly, the implication is that, while the 
anvil formed from convective detrainment, the TTL cirrus formed through some small-scale 
mechanism that involved mixing of the TTL with deep convection. 

A second example, showing possible interactions between deep convection and the TTL, 
is shown in Figs. A4 and A5 from quasi-level straight northbound flight off the coast of 
1 50 Honduras on 9 July 2002 at an altitude between 15 and 15.3 km. This case is of particular 
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interest because it shows a situation where deep convective cloud had punctured much more 
tenuous cirrus at the tropopause. The deep convective cloud sampled at 67350 and 67800 
s UTC had S — HDO values close to zero, indicative of air brought from the surface, with 
ice mixing ratio values approaching several hundreds ppmv; surrounding TTL cirrus was 
1 55 characterized by 5- HDO values of ~ —50% and ice mixing ratio values typically less than 
20 ppmv. As noted by the back-seater on the aircraft at the time, the TTL cirrus was at 
a higher level than the anvil cirrus shield that is visible in the Geostationary Operational 
Environmental Satellite (GOES)-8 imagery in Fig. A4. 

We now explore several possibilities for the origins of the TTL cirrus shown in Figure 
160 A5. 

3.1 Pre-existence 

First, the TTL cirrus may have formed prior to the convective event from some unrelated 
mechanism. Analysis of the GOES-8 imagery, matched to the four nearest pixels to the 
flight track (Minnis et al., 1998), indicates the breadth of the cloud layer visible from space 
1 65 had visible optical depths (r) sufficiently high for the cloud to act as a thermal blackbody 
(i.e. r > 10). The derived atmospheric effective temperature T, jj was ~-60 ° C where the 
WB-57F was in TTL cirrus, and ~-70 °C where it was in overshooting deep convection. 
Outside the anvil cirrus boundaries, (denoted by A and B in Figs. A4 and A5), the TTL 
cirrus also terminated, and T e j j was much higher, indicative of the lower atmosphere or 
170 surface. For comparison, in-situ MMS measurements of temperature along the extent of 
the flight leg varied about -75 °C. We interpret this discrepancy between in-situ and space- 
derived temperatures as a) indicating the TTL cirrus layer was transparent to the satellite, 
which saw only the emission from the blackbody anvil layer below, but that b) the TTL 
cirrus and anvil layers had nearly the same horizontal extent, in which case the TTL and 
1 75 anvil cirrus origins could both be attributed to the deep convection. Thus, it is unlikely the 
TTL cirrus existed previously. 

3.2 Detrainment 

If related to the convection, the TTL cirrus may simply be detrained cloudy air from the 
deep convection, that became dessicated through precipitation, leaving a tenuous layer de- 
1 80 pleted in HDO. This is unlikely. To create the abrupt gradients in 6 — HDO and Q observed 
between the convection and TTL cirrus would have required rapid precipitation of most of 
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the ice mass. Ice crystals in the convective plume were very small. Values of the ice crystal 
effective radius r e in the deep convection were derived using two techniques, from size dis- 
tributions from the CAPS probe, and from combined bulk measurements of ice water and 
1 85 optical extinction from the Harvard water and CIN probe. Bulk probes measurements of r e 
are less sensitive to condensate than the CAPS probe ( hence the incomplete data set in the 
more tenuous TTL cloud), but ostensibly they are more accurate as they do not depend on 
an assumed particle shape (Garrett et al., 2003). Both sets of measurements showed values 
of r e of about 5 to 6 pm in the the deep convective cloud. The contribution to total mass 
1 90 (measured with the Harvard probe) from precipitation ice crystals with area-equivalent radii 
larger than 25 /im (measured with the CAPS probe) was less than 10% and 15%, in the first 
and second convective penetrations, respectively. The effective radius of these precipitation 
particles, derived from CAPS measurements, was about 30 pm, which corresponds to a 
gravitational settling rate of just 10 cm s -1 . Therefore, most of the convective mass was 
1 95 in suspension; the mass that wasn’t settled slowly. Detrainment from convection could not 
have accounted for the observed TTL cirrus properties. 

3.3 Wind shear 

Alternatively, strong wind shear over the top of a mid-latitude convective domes can favor 
’’jumping cirrus” (Fujita, 1982), or breaking gravity waves that entrain ’’gulps” of convec- 
200 tive air which then spread out as thin layers of cirrus (Wang and Key, 2003). Yamamoto 
et al. (2003) showed that, in a 1 km layer above the tropical temperature minimum, Ri is 
usually less than 0.5, and sometimes less than the ’’critical” value of 0.25 below which tur- 
bulence is sustained, but that elsewhere the TTL is very stable. MMS profiles before and 
after the TTL cloud transect on July 9 indicated the cloud was embedded in a nearly isother- 
205 mal environment. Combined with measurements of wind shear just prior to and following 
the TTL cirrus flight leg, we calculated that Ri ranged from ~1.5 to ~11. Also, MMS 
measurements of the kinetic eddy dissipation rate e in the TTL cirrus were just 10 -6 W/kg. 
While this was ten times higher than what was observed in surrounding clear air, it is still 
exceptionally low. It is unlikely that vigorous cirrus formation mechanisms were playing a 
210 role in this case. 
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3.4 Gravity waves 


It has been noted previously (Potter and Holton, 1995) that deep convection can force grav- 
ity waves that propagate vertically into the stratosphere where they may induce cloud forma- 
tion. Lane et al. (2001) used numerical and analytical arguments to show that gravity-waves 
21 5 can be produced by the vertical oscillation of convective turrets about their level of neutral 
buoyancy (LNB) as they move with the wind. 

There is some evidence for convectively-forced gravity wave motions in the 9 July case 
(Fig. A5). Fluctuations in zonal wind (U) and T were found in WB-57F MMS data between 
67700 and 68300 s UTC within a layer between 15.0 and 15.3 km. The initial encounter 
220 with a convective plume at 67800 s UTC was associated with a temperature drop of ~3 K. 
As convective air was colder than its surroundings, the convective plume had risen above its 
LNB, locally squeezing isentropes upward. Two subsequent temperature drops ranged from 
1.5 to 2.2 K, but were not within dense cloud. As such, these latter fluctuations represent 
gravity wave motions forced by the convection, with the strongest within the decelerating 
225 plume itself. Nearly constant temperature during climbs before and after the TTL cirrus 
measurements indicate the environment was locally quasi-isothermal. The tropopause alti- 
tude (14.9 km) was just slightly below the aircraft track. Because gravity waves were forced 
within relatively stable air overlying the decelerating updraft, we conclude the WB-57F was 
sampling near the forcing level for the wave field. 

230 MMS time series of U and 6 perturbations show that gravity waves were first observed 
at 67000 s UTC, prior to entering the first convective plume at ~67300 s UTC, and again 
at 67800 s UTC, when the second updraft core was reached. A similar signature is seen 
in the saturation and water vapor mixing ratios. During the first encounter, only modest 
cooling was associated with a portion of the plume itself. Much of the wave energy appeal's 
235 to have dispersed in the direction of the approaching aircraft, suggesting the plume had 
significantly decayed, in contrast, cooling within the second updraft was much stronger, 
and better correlated with the total water and b-HDO time series. Since this cooling took 
place at the leading edge of the gravity wave field, we believe this particular updraft was 
forcing the observed waves that followed its penetration. 

240 At 67930 s UTC, 100 s period (in flight time) oscillations in MMS measurements of U 
and 6 appeal' to be almost in phase. This suggests waves were propagating up and away from 
the convection, but at a frequency more typically associated with decaying waves (Dean- 



Day et al., 1998). Because a low-stability region was evident in the MMS temperature 
profile above the convection at 16 km altitude, wave reflection may have also altered the 
245 phase relationship between U and 9 from the canonical quarter-cycle often observed at 
lower intrinsic frequencies. 

Because of the damping effect of the aircraft mass, MMS vertical winds are most re- 
sponsive at shorter observed periods (<50 s). Since we arc investigating somewhat longer 
periods, we instead estimate the updraft velocity of the wave from the sampling geometry, 
250 estimates of convection velocity, and linear wave theory. These result in a gravity wave 
of ~7 km wavelength having an updraft velocity amplitude of W ~ 4 m s _1 , and an in- 
trinsic frequency of Co = 0.015 s -1 (see Appendix 1). These results arc consistent with 
others outlined by Garrett et al. (2004), which described the dynamics of a TTL cirrus as 
being driven by high frequency waves with an amplitude and wavelength consistent with 
255 forcing from uplift of the TTL by a turret of deep convection. Convection produces dis- 
persive waves covering a wide spectrum in frequency and energy. However, theoretically, 
’’stationary-phase” solutions to an initial value problem for the temporary displacement of 
stratified air (Lighthill, 2001) indicate the highest frequency waves, approaching the buoy- 
ancy frequency N, have the largest vertical amplitude. Therefore, it is these that are most 
260 likely to produce the cooling required for cloud formation. 

However, while it is quite plausible that these rapid, high-amplitude gravity-wave mo- 
tions could have forced cloud formation within clear TTL air, the magnitude of the temper- 
ature fluctuations was nonetheless too small to account for the elevated levels of Q observed. 
The TTL cirrus was enriched in water by several tens of ppmv relative to surrounding air 
265 while adiabatic vertical displacement of a parcel at this level corresponds to apparent en- 
richment of just 1.6 ppmv K -1 . Second, the TTL cirrus exhibited no apparent correlation 
between the magnitude of Q, and fluctuations in 9 or U, which suggests the two were largely 
unrelated. While gravity waves may have played a role in the formation of ice mass, it was 
apparently minor. 

270 3.5 Mixing with clear TTL air 

A more likely explanation for the elevated levels of Q observed is that successive pulses of 
deep convective cloudy air mixed ice and water into the surrounding TTL. If the surrounding 
air were sub-saturated, as it might be at lower altitudes, the mixing would leave a halo of 
high humidity around the cloud. On 9 July, however, the TTL was initially supersaturated 
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275 with respect to ice, with average values of Si of about 1.2, so ice did not evaporate. Also, 
because the ice crystals were generally small, they did not settle either, and lingered as 
tenuous cloud. 

To test this hypothesis, we create a simple mixing model based on the assumption that 
total water Q, temperature T, and ice crystal number N mix linearly between convective 
280 and cloudy air, i.e., 


C mixture — (1 /) Cclear T f Cconvection 


where, £ is the mixed scalar quantity and / the fraction mixed in from convection. The mix- 
ing ratio of ice % is determined from the difference between Q mixture and SiW S at{T m i x ture), 
where w sat is the saturation mixing ratio. Cloud particle effective radius r e is proportional 
285 to the cube root of x/N- 

Fig. A6 shows the expected size and ice mixing ratio for mixed particles as a function 
of / and Si in TTL air. For the observed value of Si = 1.2, \ decreases in proportion to 
the proportion of clear TTL air mixed into the convection. However, in an environment of 
enhanced available condensate, the mixed particles grow slightly. For comparison, there 
290 is no change to modeled r e for Si = 1, and for sub-saturated air, the particles shrink then 
evaporate as / becomes small, leaving a ’’halo” of elevated humidity around the convective 
core. 

These results suggest an apparent inconsistency with the observations shown in Fig. A5: 
where measurements of x in TTL cirrus would suggest / was less than 0.1, the correspond- 
295 ing modeled value of r e in TTL cirrus would be expected to be about 0.5 gm higher than in 
the convective plume. Two sets of observations show that, in fact, r e was smaller, and by 
as much as a factor two, ranging from 2 to 4 //m. Given the generally small size of the ice 
crystals, this discrepancy can not easily be explained by appealing to precipitation removal: 
as mentioned previously, the particles were too small to have settled out rapidly. 

300 Mixing between deep convection and initially clear supersaturated TTL air easily ac- 
counts for observed ice water mixing ratios in the TTL cirrus. However it can not explain 
the observed ice crystal sizes. 


3.6 Mixing with TTL pileus cloud 

A final explanation considered for the observations is a modification to the aforementioned 
305 mixing and gravity wave hypotheses. While mixing occurred between the convective plume 
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and TTL air, the TTL air was not initially clear. Rather, because the TTL was initially su- 
persaturated, a pileus cloud formed easily ahead of the rising pulse of deep convection. The 
more energetic, and faster moving convection punctured the TTL cloud it formed, mixing in 
convective ice mass. Interfacial mixing occured between the deep convection and ambient 
31 0 air throughout its trajectory, but only where it passed through intially supersaturated layers 
of clear air did both pileus cloud form, and mixed convective ice crystals avoid evapora- 
tion. Through successive pulses of convection and mixing, TTL cirrus was formed with the 
microphysical and gravity wave characteristics shown in Fig. A5. 

To test this hypothesis requires simulation of pileus dynamics and microphysics. We 
315 arc aware of no detailed description of the dynamics of pileus formation. However, as a 
guide, dimensional analysis suggests that, from a balance between kinetic and potential 
energy, stratified isentropes arc squeezed upward by Sz ~ W/N, where N is the buoyancy 
frequency of the TTL and W the vertical velocity of the convection. Thus, deep convection 
adiabatically cools the air it displaces by 
qW 

320 ST ~ (1) 

c p N 

where g is gravity, and c p the heat capacity of dry air. Values of W in deep convection 
commonly reach ~10 m s -1 . However, above the base of the TTL, overshooting turrets 
decelerate, so W should be somewhat smaller. To first order, it increases linearly with 
height. As a guide, we assume the described gravity waves around the convection, with a 
325 vertical velocity amplitude of 4 m s -1 , and a frequency of 0.015 s -1 , were representative 
of the initial forcing of the pileus cloud. 

An estimate of the cooling forced by the displacement is of order several degrees (Eq. 
1). Condensate forms if this cooling is sufficient to raise the saturation ratio with respect to 
ice Si to approximately 1.6 (Koop et al., 2000; Baker and Baker, 2004). Above this level, 
330 the amount of condensate that forms is determined by the Clausius-Clapeyron relation. Our 
estimates of the controlling dynamics arc used to force a microphysical parcel model (Ap- 
pendix 2), which includes homogeneous freezing of solution aerosol, ice crystal diffusional 
growth, and mixing. Based on the microphysical measurements in the TTL cirrus shown in 
Fig. A5, the parcel model is intialized to values of Si and T of 1.2 and 198 K, respectively. 
335 Convective air is assumed to be 3 K colder than its surroundings, with an ice mixing ratio 
of 150 ppmv and values for r e of 5 /mi. 

Two cloud formation scenarios arc considered. In the first, a pileus cloud forms due to 
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convective uplift. The pileus cloud stays separate from the convection and does not mix (i.e. 
/ = 0). In the second scenario, after the pileus cloud forms, it is punctured by water-laden 
340 deep convection, and the two clouds mix. This is parameterized with a value of / = 0.1 
to be approximately consistent with the apparent enrichment of TTL air shown in Fig. A5. 
In both simulations (Fig. A7), cloud forms when Si reaches about 1.65, following which 
it responds to gravity wave temperature oscillations, and water vapor diffusional growth to 
nucleated ice crystals. In the warm phase of the simulated gravity wave, both the mixed and 
345 non-mixed parcels become sub-saturated. However, the subsaturation is least in the parcel 
that had been exposed to mixing with deep convection. Through repeated wave cycles, the 
cloud formation process is repeated in the non-mixed parcel, but in the mixed parcel pre- 
existing condensate inhibits the high values of Si required for new nucleation. In reality, 
new cloud formation would probably be restricted to the initial impulse in both cases, due 
350 to the dispersive nature of the wave. 

Fig. A8 shows the simulated microphysics for the mixed and non-mixed scenarios. Be- 
cause the cooling was rate was exceptionally rapid, all available aerosol in the model froze, 
yielding an ice crystal concentration of 1263 cm -3 . Karcher and Lohmann (2002) describe 
similarly high concentrations for the vertical velocities we modeled. 5 ppmv of ice water 
355 condensed out at the apex of the wave and the simulated pileus ice crystals grew to just ~ 
0.5 /jm radius. These ice crystals sizes are surprisingly small, but reflect the assumed TTL 
aerosol concentration and the high calculated value for c o; if the aerosol concentration were 
lower, the ice crystals would be larger in proportion to the cube root of the concentration 
difference; added simulation showed that for r e to be closer to values measured in the deep 
360 convection, Co would need to be one hundred times slower. 

In the non-mixed scenario (Figs. A8a and A8c), where the air parcel enters the warm 
phase of the gravity wave, it becomes sub-saturated, and the pileus ice crystals sublimate 
entirely. The cloud formation process is repeated in wave subsequent cycles. In the mixed 
scenario, where convective cloud is mixed in to the pileus layer (Figs. A8b and A8d), the 
365 pileus layer does not entirely evaporate. This is because, as deep convection mixes with its 
surroundings, it contributes not just the water vapor contained in cold dry air (as suggested 
by Sherwood and Dessler (2000)) but water from small ice crystals as well. Because these 
particles sublimate along with the original pileus ice crystals in the sub-saturated portion of 
the wave, they provide a reservoir of water vapor that limits the full extent of sublimation 
370 in the newly formed wave cloud. The cloudy air mass that is derived from TTL air is only 
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ever partially sublimated (Fig. A8d), and the cloud is sustained over repeated wave cycles. 

In Fig. A8d, the combined effective radius of the modeled pileus and convection ice 
crystals was 2 jim, which is in reasonably good agreement with the CAPS observations in 
the TTL cirrus shown in Fig. A5. While still speculative, of the various hypotheses dis- 
375 cussed for TTL cirrus formation on July 9, it appeal's that the scenario of deep convective 
air mixing with pileus cloud is most consistent with the measurements. It accounts for both 
the observed TTL condensate mixing ratios y, and the observation that the TTL ice crystal 
effective radius was smaller than in nearby convective air. In essence, neither mixing nor 
rapid gravity wave motion alone explain the measurements; it is the combination that is 
380 required. What the parcel model does not reproduce well, however, is the observed equilib- 
rium values for S t of approximately 1.2 (Fig. A7). Why S r should be greater than unity is 
unknown, and it suggests that some important physics is likely absent from our simulations. 
If our simple mixing model is forced to account for values of S, = 1.2 in the deep convec- 
tion and surrounding TTL cirrus, the results are similar (Fig. A9). If / < 0.1, as suggested 
385 from measurements of y in the TTL cirrus, the corresponding values of r e are smaller than 
2.5 pm. 

3.7 Conceptual model 

Based on these analyses, we show in Fig. A10 a schematic of three possible scenarios 
by which we suspect pileus clouds might form and evolve. In the first, (Fig. AlOa) an 
390 initially flat isentropic surface has a value of S, greater than unity ( S sa t ) but below its 
homogeneous freezing point Sf, were it to be lifted adiabatically. Vertical perturbation 
of the surface by rising deep convection causes a cloud to form above the point where it 
reaches Sf. Once the deep convection subsides to its LNB, the isentropic surfaces are left 
to oscillate about their initial position. In the warm phase of the wave, the cloud evaporates 
395 where it dips below S sa t. The cloud does not reform because the wave energy propagates, 
diminishing the wave amplitude in subsequent cycles. A second possibility (Fig. AlOb), 
however, is that the newly formed pileus cloud never evaporates. While the wave warms, it 
never warms sufficiently to dip below S sa t, and hence evaporate the cloud. In this case, the 
newly formed cirrus will be enriched in water relative to its environment only in proportion 
400 to displacement from its equilibrium level (about several ppmv per degree depending on 
the temperature). Third (Fig. AlOc), once the pileus cloud forms, it is punctured by the 
convective turret that formed it. The turret, being turbulent, mixes with the more laminar 
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pileus cloud. Ice crystals arc exchanged at their interface. Sufficient moisture is donated 
to the pileus cloud to enable it to survive the warm phase of the ensuing oscillations, even 
405 where it dips below the level that would normally be associated with S sa t in the absence of 
mixing. It is this scenario that we believe best explains the physical characteristic of the TTL 
cirrus shown in Fig. A5. Not shown is the possibility of pileus formation from isotherms 
that arc initially sub-saturated. This would require considerable lifting and cooling. Also, 
the cloud would more likely be ephemeral, as it would almost certainly evaporate in the 
41 0 warm phase of the wave. Note, also, that this diagram illustrates only individual pulses of 
convective thermals and pileus, and not rather a continued sequence of impulses that would 
blend together to form a more continuous sheet of TTL cirrus. An analogous situation is 
one where a series of convective bubbles reaching the tropopause relax to their LNB to form 
continuous laminar anvil cloud. 

415 4 Susceptibility for pileus formation in the tropical TTL 

While it is evident from Fig. A1 that pileus clouds can form in the TTL, it remains un- 
known how often the TTL might be susceptible to pileus formation above convection. The 
likelihood of pileus formation at any location depends on three factors: the stratification of 
the TTL, the velocity of the convective uplift, and foremost, the humidity: higher values of 
420 Si require less lifting by convection to initiate ice crystal nucleation. In the tropics, at least, 
regions of TTL cirrus arc highly correlated with Si, and arc most common in regions where 
deep convection is active (Sandor et ah, 2000; Dessler and Yang, 2003; Wu et al., 2005). 
We have no evidence to show that these TTL cirrus layers might originate as pileus cloud. 
However, it should be considered a possibility. 

425 To illustrate. Fig. All shows Microwave Limb Sounder (MLS) measurements obtained 
during 1992 aboard the Upper Atmosphere Research Satellite (UARS) (Read et ah, 2004) 
The frequency curves show that air is often supersaturated in the lower tropical TTL, but 
only rarely so aloft (Fig. All). Superimposed on these curves is the amount of isen- 
tropic lifting that would be required to cool air to the point where ice nucleation would be 
430 expected, in accordance with the Koop et al. (2000) formulation. For air that is initially 
saturated with respect to ice, the required lifting is about 300 m, independent of altitude. 
To produce such uplift above deep convection, assuming a balance between kinetic and po- 
tential energy, W ~ N6z, the associated deep convective updraft velocity would be 4 m 


14 



s _1 . This seems reasonable given that in the TTL, convection is above its LNB, and it is 
435 therefore decelerating. Assuming such lifting, then at 146 hPa, in the lower TTL, ~ 20% 
(10%) of land (ocean) air would be expected to form pileus cloud over deep convection. 

This estimate of susceptibility of the TTL to pileus formation must be considered with 
caution given that the uncertainty in MLS calculations of S) is approximately 30%, 40% and 
50% at 146 hPa, 121 hPa, and 100 hPa, respectively. It is probably conservative, because 
440 MLS values of .S', represent an average over a vertical depth of ~4 km, whereas high hu- 
midity layers arc often thinner (and presumably more common). Higher resolution aircraft 
measurements off Honduras and Costa Rica show air mostly supersaturated into the middle 
TTL (Fig. A1 1). Si mi lar results have been noted in TTL air near Costa Rica (Jensen et ah, 
2005). 

445 5 Conclusions 

The relevance of the formation of pileus cloud to how phases of water arc partitioned in the 
TTL depends in large paid on their longevity. Ordinarily it would be expected that the cloud 
would exist only in the phase of the gravity wave where S', > 1. Once the potential energy 
associated with the vertical displacement of the initial impulse has dispersed, the cloud dis- 
450 appeal's. What we have argued here is that, under some situations, pileus will linger past 
the duration of the convective impulse. Overshooting deep convective clouds may be dry 
compared to their warmer surroundings, but they are also laden with water in the form of 
small ice crystals: their tops would not be clearly visible otherwise. Because convection 
is turbulent, it necessarily mixes into its surroundings both the ice and the water vapor. If 
455 the environment is sub-saturated, the ice evaporates where we normally interpret the edge 
of the cloud to exist, leaving behind a halo of humid air. If the environment is supersatu- 
rated, however, it is probable that uplift ahead of the convection will be sufficient to initiate 
condensation. Where this pileus cloud is punctured by the faster moving convective turret, 
ice crystals are mixed from the convection into a cloudy atmosphere and do not completely 
460 evaporate. In fact, as described by our parcel model simulations, these ice crystals act as 
a reservoir of water vapor, which inhibits deep subsaturation in the warm phase of ensuing 
buoyancy oscillations. As a result, the original pileus condensate never fully evaporates, 
even once the convection has subsided. A repartion from water vapor to ice has occured 
within the fraction of the mixed air that was derived from the previously clear TTL. 
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465 In examination of the origins of a TTL cirrus layer formed off the coast of Honduras, it 
was this mechanism, of several that were considered, that appeared most consistent with 
the observations. The ultimate fate of the TTL cirrus ice crystals in this case is unknown. 
However, other observations of TTL cirrus associated with deep convection suggest the 
clouds to be relatively long-lived. For example, a study of the evolution of a TTL cirrus 
470 layer above an anvil during CRYSTAL-FACE showed that the cloud properties were stable, 
even while the anvil beneath dissipated (Garrett et al., 2005). The longevity may have been 
partly because the ice crystals were small and in a supersaturated environment, so that they 
were unlikely to precipitate or evaporate. The TTL cirrus ice crystals were only several 
micrometers radius, in which case gravitational settling was only several tens of meters 
475 per day. Also, it has been suggested Hartmann et al. (2001), and demonstrated by Garrett 
et al. (2005), that TTL cirrus cools when a cold anvil layer is present below. However, such 
cooling does not appeal - to be necessary for sustenance. Comstock et al. (2002) noted TTL 
cirrus in the Tropical Western Pacific that persisted for days following a convective event, 
even once the anvil had dissipated. 

480 The results described here suggest a hypothesis in which, by forming pileus cloud in cold, 
humid environments, convection plays a lasting role in how water and ice are partitioned. 
The redistribution is due to small-scale laminar and turbulent processes at altitudes well 
above the convective level of neutral buoyancy where anvil cirrus is normally found. While 
we provide photographic evidence that pileus clouds form at high altitudes and become 
485 TTL cirrus, we know of no in-situ measurements that are conclusively within the pileus 
clouds themselves. Hence, while plausible, the hypothesis remains speculative. Space- 
borne measurements of the TTL with lidar and radar should help show where and how 
often thin cirrus layers are present above deep convection. Detailed cloud resolving model 
simulations might show how these clouds evolve, and quantify their potential influence to 
490 global fields of vapor and radiation near the tropopause. 

Al Wave sampling 

As viewed in the horizontal plane, gravity waves generated by a convective updraft will form 
a steady-state, parabolic-shaped field of constant phase curves, upstream from the energy 
source (Fig. A12). In time, each of the parabolic curves tracks the horizontal movement of 
495 the underlying convection, orienting its axis along the horizontal wind vector at the altitude 
of wave forcing, here, the wave vector k is defined to lie parallel to the mean wind, but by 
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convention must have a positive eastw ard component. 

In general, the complex geometry of a gravity wave field can only be fully determined by 
means of three-dimensional, non-hydrostatic, time-dependent mesoscale models. However, 
500 one can approximately describe the pattern formed during an in-situ encounter with a wave 
field within a limited geographical area. If the observed wave periods are relatively constant 
while sampling along a quasi-level altitude, local parabolic curvature may be neglected, and 
lines of constant phase may be assumed. 

Fig. A13 depicts horizontal sampling of a single, monochromatic gravity wave forced by 
505 a convective updraft and moving at the velocity of the underlying air column. The sampling 
aircraft encounters a wave with phase Vfi at point A. After a time interval St, the next wave 
is encountered at the same phase (i.e., ^2 = Vfi + 27r) at point B. 

A system of equations is used to determine wavelength A and phase speed c from the 
apparent wavelength A m , and period St derived from aircraft measurements: 

510 A m = V A St 

c=V c cos 0 

A = I A m COS (j) — c5t\ 

where, 6 and <p arc the angles between the wave vector and either the convection velocity 
V c or aircraft velocity Va, respectively. For the wave encounter of 9 July 2002 (Fig. A5), 
515 we assume the wave is forced by convection at or just below the flight altitude. Here, the 
wave vector k is aligned with the mean wind (U = 11 m s 1 , V = 0 m s -1 ) at the WB-57 
altitude. In this circumstance, intrinsic frequency Co is defined as, 

*= 2 i(c-U) 

On this day, the NASA ER-2 flew in stacked formation with the WB-57F over the same 
520 flight track, releasing a dropwindsonde approximately 250 km south-southeast of the con- 
vective cell, and less than 30 min prior to sampling of the wave by the WB-57F. Assuming 
the wave field travels with the convection that forces it, we present solutions for A, c, and Co 
based on the convection velocity V c from the sonde wind data. 

The first method uses the pressure-weighted mean wind, computed from the ocean sur- 
525 face (1006 hPa) up to the tropopause at ~132 hPa. Calculations yield V c = 6.1 m s~ 1 from 
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104°. Here, the wave and convection velocity vectors arc separated by 9 = 166°, resulting 
in c = —5.9 m s' 1 . From this, A = 7.2 km and Co = 0.0148 s' 1 . 

Alternately, the steering wind at 700 hPa may be used. Here, V c = 6.4 m s' 1 from 88°. 
Since the wave vector is the same, it follows that 6 = 182°, yielding c = —6.4 m s' 1 , 
530 A = 7.2 km, and Co = 0.0153 s' 1 . 

Local static conditions at flight level yield a buoyancy frequency N = 0.0164 s' 1 , so 
the waves are close to the criterion for decay {Co ~ 0.9A’). We use a mean phase speed 
(c = —6.2 m s' 1 ) and wavelength (A = 7.2 km) to determine the updraft velocity. Air 
flowing through the wave requires a transit time of 209 sec (~3.5 min) to move from the 
535 warm phase to the cold phase. Assuming convection initially displaced the parcel upward 
by 300 m, the average updraft velocity is W = 2.9 m s' 1 , with a peak velocity W max = 4.1 
ms -1 . If our calculations were in error such that the waves were formally evanescent (i.e., 
Co = N ), a slightly shorter wavelength (A = 6.6 km) would yield the strongest vertical 
motions possible {W = 3.1 m s' 1 and W max = 4.4 m s' 1 , respectively). 


540 A2 Homogeneous nucleation model 


Our model simulates the evolution of individual sulfuric acid solution particles within a 
parcel of clear TTL air. Consistent with measurements in the TTL (Lee et al., 2003), the size 
distribution of dry particles is log-normal and bimodal, with mode radii, spectral widths, 
and concentrations of 0.015 /mi and 0.06 /im; 1.4 and 1.5; and 2500 cm' 3 and 150 cm' 3 , 
545 respectively. The parcel oscillates vertically with a frequency Co and a vertical velocity 
amplitude W. Being a parcel model, environmental wind shear is ignored. 

The model permits mixing at the interface between pileus and deep convective cloud. At 
the apex of the first cycle in the wave, it is assumed this mixing occurs instantaneously. 
Water and sensible heat are mixed linearly: i.e., 


550 C mixture — (1 f) Cpileus T" f Qconvection 


where, Q is the mixed quantity and / the fraction mixed in from the convection. The model 
then tracks the sizes of individual particles as they adjust to instantaneous (in the case of 
haze) and diffusive (in the case of ice) equilibration with the rapid temperature fluctuations 
induced by convection. 

555 At equilibrium, solution aerosol adjust their volume V such that the activity a w of the 
curved solution droplet is equivalent to the relative humidity in surrounding air. It has been 
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argued by Karcher and Koop (2005) that aerosol composition effects on surface kinetics 
can retard this equilibration under conditions of rapid cooling. We believe, however, that 
quasi-equilibrium is maintained for the following reason. The time scale for diffusion of 
560 vapor to an aerosol particle uujj is of order r 2 /D, where r is the aerosol radius and D the 
diffusivity. This time scale is easily modified to account for kinetic limitations associated 
with molecular diffusion to small droplets, and by a resistance factor - the condensation 
coefficient - which accounts for the efficiency of water molecule uptake. Even if the con- 
densation coefficient is very low, t ( h f j is many orders of magnitude less than the time scale 
565 associated with changes in ambient humidity. For any plausible scenario of atmospheric 
cooling, an aerosol particle must always be in an effective equilibrium with its environment 
(Pruppacher and Klett, 1997). 

The nucleation rate of ice embryos in solution J (a w . T) determines the probability P = 
1 — exp (— JV At) an aerosol freezes in time period At (Koop et al., 2000; Baker and Baker, 
570 2004). According to this formulation, aerosol freeze when the saturation ratio with respect 
to ice S{ reaches ~ 1.6. Large aerosol have the highest probability of freezing (V ~ r 3 ), but 
smaller aerosol tend to have correspondingly high concentrations (N a ~ r -3 ). Therefore, 
it is almost equally likely that aerosol of a given size will homogeneously freeze during At. 

Once frozen, aerosol grow rapidly by vapor diffusion. This depletes the ambient vapor 
575 field and therefore a w \ further nucleation may be suppressed so that only a fraction of 
available aerosol become ice crystals. The assigned value of the condensation coefficient 
is 0.2 (Delval and Rossi, 2004). High ice crystal concentrations are favored by low T and 
high W, and are only weakly related to the size distributions of the aerosol (Karcher and 
Lohmann, 2002). 
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Figure Captions 


Fig. Al. Photo of pileus forming on top of deep convection over land near Darwin, Australia on 14 
Nov 2004 at 7:05 pm local time. 


Fig. A2. Photographs of pileus and thin cirrus formation obtained at 12 km altitude over Louisiana 
on 30 April 2004 during the Mid-Latitude Cirrus Experiment (MidCiX). Photograph (B) is a closeup 
of the convective dome shown 6 min earlier in (A). 


Fig. A3. Profile of (A) temperature T and potential temperature 9 , (B) H 2 0 mixing ratios (measured 
vapour: black; saturation vapour with respect to ice: grey; total (vapour plus ice): red) and J-HDO 
in a profile obtained on 28 July, 2002 over southern Florida. Two isotopically distinct cloud layers 
are seen: anvil cirrus outflow from convection located below 13 km altitude, and a more tenuous 
TTL layer centered at 14.5 km altitude 


Fig. A4. WB-57F flight leg through cloudy air at approximately 15 km altitude off the East coast of 
Honduras during the CRYSTAL-FACE experiment, in 5 minute segments beginning at A on 66300 
s UTC 9 July 2002. Arrows point to regions of overshooting deep convection. 

Fig. A5. In-situ and satellite retrievals corresponding to the flight leg shown in Fig. A4, and bounded 
by letters A and B. Measurements show, GOES-8 retrievals of cloud effective temperature T e ff and 
optical depth r, in-situ perturbations to environmental potential temperature A 6 and zonal wind 
velocity A(7 from the MMS instrument, the saturation ratio with respect to ice 5) and total Q and 
vapor uj water mixing ratios from the Harvard water probe, fractionation of HDO relative to H 2 0 
8 — HDO from the ALIAS probe, and ice crystal effective radius r e from the CAPS probe (black 
line), and the CIN and Harvard water probes (red dots). Shaded area shows mixing ratios of ice. 


Fig. A6. Particle radius (/xra) (red) and ice water mixing ratio (ppmv) (blue) for a hypothetical 
mixing scenario between the convective cloud shown in Fig. A5 and ambient TTL air, for values of 
clear air TTL S) shown. 


Fig. A7. Simulations of .S’, in a simulated non-dispersive gravity wave for two scenarios: (blue) 
pileus cloud forms, and then reforms in successive pulses of a gravity wave, and (red) the pileus 
cloud forms, then mixes with a deep convective plume. Mixing occurs at 100 s. 
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Fig. A8. As in Fig. A7 for non-mixed (A and C) and mixed (B and D) scenarios. A and B rep- 
resent simulations of w (dashed line), Q (solid line), and \ (shaded area), and C and D represent 
simulations of particle radius (shaded area indicating bounds of one standard deviation away from 
the mean). 


Fig. A9. As in Fig. A6, except convective air is mixed in fractional quantity / into pileus cloud with 
X = 5 ppmv and r e = 0.5 /im. 


Fig. A10. A schematic diagram showing three scenarios for the hypothesized formation of pileus 
cloud by convection, and the subsequent evolution of the forced gravity wave. See discussion in text. 


Fig. All. Distributions of S) in the tropical TTL retrieved 20° S and 20° N using Version 7.02 
Microwave Limb Sounder (MLS) measurements aboard the UARS satellite (lines), and in situ values 
from two flights off Honduras flown during the July 2002 CRYSTAL-FACE project (shaded area). 
Contours represent the estimated isentropic lifting (in meters) required to initiate freezing in air with 
the given pressure and ,S’ ( . 

Fig. A12. Parabolic pattern formed by gravity waves generated by a convective source. Viewed 
from above, the updraft forcing the wave field is indicated by a black dot; two-dimensional curves of 
constant phase L are shown to intersect and lie upstream from the convection. Va indicates aircraft 
velocity along the ground-based track (bold line); the wave vector k is shown as the axis of the 
parabolic phase curves at time t. 


Fig. A13. Geometry of wave sampling, as viewed from above aircraft flight track, tpi and -02 are 
phase lines of a monochromatic gravity wave sampled at times t\ (solid) and < 2 (dashed), respec- 
tively. Points A and B show the position of the aircraft during consecutive wave samples. Distance 
between phase lines equals the wavelength A, and the wave vector k is shown orthogonal to the 
phase lines. Va and V c are aircraft and convection velocities with respect to earth, o and 6 are an- 
gles between the wave vector and either the aircraft track or convection velocity vector, respectively. 
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Fig. Al. Photo of pileus forming on top of deep convection over land near Darwin, Australia on 14 
Nov 2004 at 7:05 pm local time. 


A 

TTL cirrus 

B Zoom 6 min later 

TTL cirrus 




Pileus 











Anvil 

Convection 

Anvil 



Pileus 


Fig. A2. Photographs of pileus and thin cirrus formation obtained at 12 km altitude over Louisiana 
on 30 April 2004 during the Mid-Latitude Cirrus Experiment (MidCiX). Photograph (B) is a closeup 
of the convective dome shown 6 min earlier in (A). 
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Fig. A3. Profile of (A) temperature T and potential temperature 0, (B) H 2 O mixing ratios (measured 
vapour: black; saturation vapour with respect to ice: grey; total (vapour plus ice): red) and (5-HDO 
in a profile obtained on 28 July, 2002 over southern Florida. Two isotopically distinct cloud layers 
are seen: anvil cirrus outflow from convection located below 13 km altitude, and a more tenuous 
TTL layer centered at 14.5 km altitude 
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Fig. A4. WB-57F flight leg through cloudy air at approximately 15 km altitude off the East coast of 
Honduras during the CRYSTAL-FACE experiment, in 5 minute segments beginning at A on 66300 
s UTC 9 July 2002. Arrows point to regions of overshooting deep convection. 
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Fig. A5. In-situ and satellite retrievals corresponding to the flight leg shown in Fig. A4, and bounded 
by letters A and B. Measurements show, GOES-8 retrievals of cloud effective temperature T e ff and 
optical depth t, in-situ perturbations to environmental potential temperature Ad and zonal wind 
velocity A U from the MMS instrument, the saturation ratio with respect to ice .S', and total Q and 
vapor w water mixing ratios from the Harvard water probe, fractionation of HDO relative to H 2 O 
6 — HDO from the ALIAS probe, and ice crystal effective radius r e from the CAPS probe (black 
line), and the CIN and Harvard water probes (red dots). Shaded area shows mixing ratios of ice. 
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Fig. A6. Particle radius (red) and ice water mixing ratio (ppmv) (blue) for a hypothetical 

mixing scenario between the convective cloud shown in Fig. A5 and ambient TTL air, for values of 
clear air TTL S; shown. 
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Fig. A7. Simulations of Si in a simulated non-dispersive gravity wave for two scenarios: (blue) 
pileus cloud forms, and then reforms in successive pulses of a gravity wave, and (red) the pileus 
cloud forms, then mixes with a deep convective plume. Mixing occurs at 100 s. 
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Fig. A8. As in Fig. A7 for non-mixed (A and C) and mixed (B and D) scenarios. A and B rep- 
resent simulations of w (dashed line), Q (solid line), and x (shaded area), and C and D represent 
simulations of particle radius (shaded area indicating bounds of one standard deviation away from 
the mean). 
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Fig. A9. As in Fig. A6, except convective air is mixed in fractional quantity / into pileus cloud with 
X = 5 ppmv and r e = 0.5 fim. 
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Fig. A10. A schematic diagram showing three scenarios for the hypothesized formation of pileus 
cloud by convection, and the subsequent evolution of the forced gravity wave. See discussion in text. 
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Fig. All. Distributions of Si in the tropical TTL retrieved 20° S and 20° N using Version 7.02 
Microwave Limb Sounder (MLS) measurements aboard the UARS satellite (lines), and in situ values 
from two flights off Honduras flown during the July 2002 CRYSTAL-FACE project (shaded area). 
Contours represent the estimated isentropic lifting (in meters) required to initiate freezing in air with 
the given pressure and Si. 
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Fig. A12. Parabolic pattern formed by gravity waves generated by a convective source. Viewed 
from above, the updraft forcing the wave field is indicated by a black dot; two-dimensional curves of 
constant phase '0 are shown to intersect and lie upstream from the convection. Va indicates aircraft 
velocity along the ground-based track (bold line); the wave vector k is shown as the axis of the 
parabolic phase curves at time t. 
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Fig. A13. Geometry of wave sampling, as viewed from above aircraft flight track, ipi and 02 are 
phase lines of a monochromatic gravity wave sampled at times t,-\ (solid) and t% (dashed), respec- 
tively. Points A and B show the position of the aircraft during consecutive wave samples. Distance 
between phase lines equals the wavelength A, and the wave vector k is shown orthogonal to the 
phase lines. Va and V c are aircraft and convection velocities with respect to earth. 0 and 9 are an- 
gles between the wave vector and either the aircraft track or convection velocity vector, respectively. 
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